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Abstract 

Previously it has been shown that some classes of mixing dynamical systems have limiting 
return times distributions that are almost everywhere Poissonian. Here we study the behaviour of 
return times at periodic points and show that the limiting distribution is a compound Poissonian 
distribution. We also derive error terms for the convergence to the limiting distribution. We also 
prove a very general theorem that can be used to establish compound Poisson distributions in many 
other settings. 



1 Introduction 

In 1899 Poincare showed that for a map T on some space il which has an invariant probability measure 
/x, almost every every point returns within finite time arbitrarily close. This means that for every 
(measurable) AcSi with fi(A) > the return time function ta(x) = min{/s > 1 : T k x E A} is finite 
for /x-almost every x G A. This result was quantified by Kac in 1947 for ergodic measures. His theorem 
states that J A ta{x) dfi(x) = 1, provided \i is ergodic, which implies that ta(x) is on average equal to 
1/ n{A). Since 1990 there has been a growing interest in the statistics of return times and in particular 
in the distribution of ta- Considering that it was shown in (2Ql [19] that for ergodic measures the 
limiting distribution of a sequence of (rescaled) return functions Tu„ for n — > oo can be any arbitrarily 
prescribed distribution for suitably chosen sets U n , it is necessary to assume that the return sets A are 
dynamically regular. 

For the measure of maximal entropy on a subshift of finite type, Pitskel 25J showed that the return 
times are in the limit Poisson distributed for cylinder sets A n (x) (fi(A n (x)) — > as n — > oo) where the 
set of suitable 'centres' x G form a full measure set. A similar result had independently been obtained 
by Hirata [T5] by a different method. For ^-mixing Gibbs measures Galves and Schmitt [TO] showed 
in 1997 that the first return time is in the limit exponentially distributed and that the convergence is 
at an exponential rate. Subsequent results (e.g. [2TJ [8] [17l [13l [5] [1] ) established limiting distribution 
results for first or multiple return times in various settings and sometimes with rates of convergence. 

Almost all previous results look at the distribution of return times near generic points. The notable 
exception being the paper [15] by Hirata which gives the distribution of the first return time at a periodic 
point. In the present paper we consider periodic points for sufficiently well mixing invariant measures 
and show that the limiting distribution is compound Poissonian. The compound Poisson distribution 
has previously been used in various settings including the analysis of internet traffic where the waiting 
time between packets is exponential and the size of each packet is geometrically distributed. It has 
also been used to model the survival of capitalist enterprises in the free market system |21) . The main 
technical result, Proposition [TJ provides conditions under which one obtains a compound Poissonian 
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distribution with error terms. This result in itself will be of interest to a much larger community than 
the one of dynamicists addressed in this paper. 



Let /ibea probability measure on a space fl which carries a transformation T, preserving fi, and whose 
cr-algebra is generated by the joins A k — V^=o T^^A, k = 1, 2, . . ., of a given finite measurable partition 
A. The elements of A k are called k-cylinders. We assume A is generating, i.e. the elements of .4°° are 
single points. Denote by \A the characteristic function of a (measurable) set A and define the random 
variable: 

T 

3=1 

z € fi. The value of Qa measures the number of times a given point returns to A within the time r. 
Typically the obsevation time r is chosen to be [t//i(A)} where t is a parameter value. (The rescaling 
factor l/fi(A) agrees with Kac's theorem.) For instance, if \x is the measure of maximal entropy on a 
subshift of finite type, then Pitskel [55] showed that Ca„(x) is for /i-almost every x € £1 in the limit 
n — ► oo Poisson distributed (where r„ = [t / fi(A n (x))] and A n (x) denotes the unique n-cylinder that 
contains x). In [13] we have proven a similar result for a much wider class of systems and provided 
error estimates. 

We develop a mechanism which allows to prove the compound Poisson distribution of return times 
at periodic points x and also to obtain error estimates as the cylinder sets A n {x) shrink in measure to 
zero. 

To be more precise, if we denote by Cn( z ) the counting function Y^j=iXA n (x) ° T^{z), with the 
observation time t„ = h^JJTa (x)) > we Wlu " study the following distribution: 

P(C=r), r = 0,1,2,..., (1) {DE} 

where t > is a parameter and p £ [0,1) depends on the periodic point x and will be given in Sect. 3. 
We will show that the limit n — > oo is the compound Poisson distribution (see Section 2) if /j, is a (cf>, /)- 
mixing measure. We also provide rates of convergence. This then implies under some mild additional 
assumptions [2] the uniform integrability of the process £*. 

We then extend this result to return times, i.e. to the distribution of the process Cn( z ) restricted to 

the cylinder A n (x). The measure /i is then replaced by the conditional measure fi n = ^^ /j 

We refer to this second case as the distribution of the number of visits for return times. 

Our results for return times considerably improve on the work of Hirata [15] , where he computed 
(without error) the distribution of the first return time (order r — 0) around periodic points and for 
Gibbs measures on Axiom- A systems. 

The plan of the paper is the following. The purpose of section 2 is to prove Proposition Q] that gives 
general conditions under which a sum of mutually dependent 0, 1-valued random variables converges to 
the compound Poisson distribution and provides error terms. A similar result that had been inspired 
by a theorem of Sevast'yanov |26j . was proved and used in [13j for the Poisson distribution. 

The distribution of return times is tied to the mixing properties of the invariant measure considered. 
For that purpose we introduce in the third section the (<p, /)-mixing property. This property is more 
general that the widely used ^-mixing property and is reminiscent of Philipp and Stout |24j 'retarded 
strong mixing property'. In this way one can obtain distribution results on return times of some 
well studied dynamical systems that are not 0-mixing, e.g. rational maps, parabolic maps, piecewise 
expanding maps in higher dimension .... 

The third section is devoted to the proof of the existence of the limit distribution and rates of 
convergence for entry times (Theorem 7), while the fourth section extends those results to return 
times (Theorem 10). Section 5 contains a careful application to rational maps with critical points 
(Theorem 11). 

We conclude this introduction with an interesting observation. Limit distributions for entry and 
return times have been provided along nested sequences of cylinder sets converging to points x which 
were chosen almost everywhere or as periodic points. In section 3.4 we will show how to find points 



A n {x) 
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x which do not have limit distributions at all, and this will be achieved by using our results on the 
compound Poisson distribution around periodic points. 

2 Factorial moments and mixing 

The main purpose of this section is to prove a very general result which we use to prove the main 
results in sections 3 and 5 but which can also be useful to establish compound Poisson distribution 
with respect to the geometric distributiion in many other settings. For more general compound Poisson 
distributions see [9]. More recently (e.g. [HE]) there have been efforts to approach compound Poisson 
distributions using the Chen-Stein method. The treatment in [7J has a more general setting, but the 
result is far from applicable to our situation^ Proposition [T] is of general interest and is reminiscent 
of existing theorems which establish the Poisson distribution (cf. [MJ from the convergence of 

the moments. It provides general conditions under which the distribution of a finite set of 0, 1-valued 
random variables is close to compound Poisson (and provides error terms). In sections 3 and 5 we 
then use it to obtain the speed of convergence for the limiting distributions for ^-mixing systems, some 
non-Markovian systems and equilibrium states for rational maps with critical points. 

2.1 Compound Poisson distribution 

For a parameter p € [0, 1) let us define the polynomials 

r = 1,2,..., where Pq = 1 (r = 0). The distribution e^ t P r (t,p), r = 0, 1, 2, ... is sometimes called the 
Polya-Aeppli distribution fTS]. It has the generating function 

oo 

g p (z) = e~* ^2 zTp r = e* 1 ^ , 

r=0 

its mean is and its variance is t n~j[p ■ Note that for p = we recover the Poisson terms e~ l P r (t, 0) = 

e — * -^j and the generating function go(z) = e*^ 1 ) which is analytic in the entire plane whereas for 
p > the generating function g p (z) has an essential singularity at i. The expansion at zq = 1 yields 

9 P { Z ) = J2T=a( z ~ 1 ) fe( 3fe where 

(Qo = 1) are the factorial moments. Note that in particular Po(0,p) = 1 and P r (0,p) — for all r > 1. 

2.2 Return times patterns 

{returntimespatt 

Let M and m < M be given integers (typically m « M) and let r £ N be some (large) number. For 
r = l,2,3,...we define the following: 

(I) G r (r); We denote by G r (r) the r- vectors v = (v%, . . . , v r ) £ Z r for which 1 < v\ < V2 < ■ ■ ■ < v r < r. 

(II) G r> j(r): We divide the set G r into disjoint subsets G rj where G r j consists of all v £ G r for which 
we can find j indices i\, 12, ■ ■ . , ij £ {1, 2, . . . , r}, i\ = 1, so that Vk — Wjt-i < M if k ^ ■ ■ ■ , ij and so 
that Vk — Vk-i > M for all k = 12, ■ ■ ■ , ij- 

For v £ G r j the values of Vi will be identified with returns; returns that occur within less than 
time M are called immediate returns and if the return time is > M then we call it a long return (i.e. 
if Vi+ 1 — Vi < M then we say v i+ i is an immediate return and if Vi + \ — V4 > M the we call V4 a long 
return) . That means that G r _j consists of all return time patterns v which have r — j immediate returns 

x We thank the referee for pointing us towards Chen and Roos' work and also for other enlightening remarks. 
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that are clustered into j blocks of immediate returns and j — I long returns between those blocks. The 
entries v ik , k = 1, . . . , j, are the beginnings (heads) of the blocks (of immediate returns). We assume 
from now on that all short returns are multiples of m. (This reflects the periodic structure around 
periodic points, cf. condition (II) of Proposition [TJ) 

(III) G r ,j iW (r): For v G G r .j the length of each block is v ik+1 _i — Vi k , k = 1, . . . ,j — 1. Consequently 
let us put Wk — 7^(vfc — ffc-i) for the individual overlaps, for k ^ i\, i 2 , ■ ■ ■ , ij- Then Y^ k J^ +1 we = 
^( v i k +i-i ~ v ik) 1S the overlap of the fcth block and w = w(v) — X^fc^ii i 2 i Wk ^ ne overlap of 
v. We now put G r .j iW ={»£ G r j : w(v) = w}. (G r j = \J W G r .j, w is a disjoint union.) 

(IV) A(v): For v in G r j we put 

A(v) = mm{v lk - v ik _ x : k = 2,...,j} 

for the minimal distance between the 'tail' and the 'head' of successive blocks of immediate returns (or 
the length of the shortest one of the long gaps) . 

2.3 Compound Poisson approximations 

The purpose of this section is to prove the following result on the approximation of the compound 
Poisson distribution. 

{sevastyanov} 

Proposition 1 Let M,m,T be as above. Let rjj, j = 1, ...,r, be 0,1-valued random variables on 
some fl for v G G r put 7^7 = rj Vi . Choose 5 > and define the 'rare set' R r = \Jjj=iRr,j; where 
Rr.j = {v G G r .j : A(v) < 5}. Let /x be a probability measure on O which satisfies the following 
conditions (Co is a constant): 

(L) E(r)j) = (3 for all j = 1, . . . ,r (invariance of the measure). 

(LL) Suppose that there are numbers < p_ <p< p+, 4> > so that for all v G G r .j iW \ R r .j 

|Efar) - P w P j I < Q> W - P-) + P w ((1 + - !) 
if all of the individual overlaps wi are multiples of m, and 

E(w) = 

if some of the individual overlaps we are not multiples of m. 

(LLL) There are some constants 7 > 1,71,72 small (e.g. 7(71 + 72) < j^), so that for all r 

S E ( , s ,< c ^ti:(i:!) 7 r^(;::!)7r». 

veR r J=2 s=l v ' 

Let us put ( = Vj an d t = (1 — p)t(3. 

Then there exists a constant C\ so that for every t > one has 

|P(C - r) - e-*P r (t )P )\ < C 1 ( 11 +5[3)t r - 1 — + C 1 (p? + P+ - P -+q>) 

Note that the constants 7 and 72 don't enter the final estimate in an explicit way. The significant 
quantity here is 71 which typically is « 1 where 7,72 only have to be small enough of order 0(1). 

The choice of S is central to the application of this proposition. In the application however <f> depends 
on S and in fact cf>(5) — > 0+ as S — * 00. Obviously a larger value for S increases the error term as one 
sees from the expression, but also a smaller value increases the error term since the rare set R r becomes 
larger and the 'mixing property' in (II) will require larger p + and smaller p_, thus again increasing the 
error estimate. The trick is to optimise 8. 

Proof. We compare the generating function ip(z) for the process ( with the generating function g p (z) 
for the compound Poissonian. In part (A) we compare their Taylor coefficients at z = 1 and in part 



te 2r+ 



t< - if 

, r t 1 + 2p -r 

(2p) r e 1 - i p if 



t > 
t < 



\pr 
\pr 
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(B) we use Cauchy estimates to compare their Taylor coefficients at z = which then gives us the final 
result. 

(A) The coefficients at z = 1 (factorial moments) of the generating function tjj(z) = ^^{C = 

r ) = E^Lo( z - 1 Y lJ r are 

U r = E (^)> 

tfeG r 

while the coefficients of the generating function g p (z) = Y^Lo( z ~^) r Qr{t,p) f° r the- compound Poisson 
distribution are 

where i = (1 — p)/?T. We will now compare the coefficients Q r to the coefficients U r . There are three 
parts to the comparison: (i) Assumption (II) is used to compare the terms for which v e G r \ R r ; (ii) 
Assumptions (I) and (III) are used to estimate the total contributions made by v G R r ; (iii) We have 
to estimate the contribution to Q r that correspond to overlaps u which do not occur for vectors v in 
G r and therefore cannot be matched with terms in the sum that defines U r . 
More precisely, we estimate as follows: 

\U r - Q r ((l-p)(3r,p)\ < EE E l E (^) -P W ?\ + E ( E (^) +P ui * ) P i ) + V(r). 

j u v£G r j yU \Rj veR r 

Before we proceed to bound the three terms on the right hand side let us estimate the cardinality of 
the sets G rj >- (Note that u > r — j if G rj > is nonempty.) Since G rjjtl consists of all v e G r j that 
have a total overlap u (in j blocks of immediate returns) we get 

ii-(r-)) + r- )- l \ / r — 1 \ _ r J / u — 1 \ f r-1 \ 
r-j-l ' )\ j-1 J \ r-j-1 J { j-1 J 

(j blocks positioned 'anywhere' on an interval of length r, u overlaps distributed on r — j immediate 
returns and j blocks beginning on any of the r return times). 

Now we estimate the three terms in the coefficient comparison as follows: 
(i) The first error term (difference between the dominating terms) is bounded using assumption (II): 

r oo 

EE E Hv*)-P w P j \ 

j=l u=r-j v£G r ,j,u\Rr 
r oo 

< E E IGwlPiPl-ptii-fl) 

j=l u=r-j 




< c 2 q{r - l)Q r -i(t,p+) + <j)Q r (t,p + ) 



(because {r — j) y ^ _ \ j = (r — 1) y r - _ \ j), where t = (1 — p)t(3 and q = p + — p_. 

(ii) For the second term let us note that R r = (J . R r j where R r j = {«£ G r j : A(v) < 5}. Put i?* 
for those v G R r where v i+1 — Vi > S for exactly s — 1 indices i\, 12, ■ ■ ■ , i s -i and put i s = v r (obviously 
1<s<j — l<r — 1 and i s _i < r — 1). 



|G r .j jtt | < 
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To estimate the cardinality of i?* ■ u — R r ' n G r ,j.u let us note that the number of possibilities of 
v il < v i2 ■ ■ ■ < v is (entrance times for long returns bigger than 5) is bounded above by ^ r s (this is the 
upper bound for the number of possibilities to obtain s — 1 intervals contained in the interval [1,t]), 
and each of the remaining j — s return times less than 6 assume no more than 5 different values. Since 
the indices i s , . . . , ik s can be picked in ^ lZ\ ) many ways out of j blocks, we obtain: 



j 



1 \ 5 j - a / r - l \ / u - 1 



W - y s _ i j s! \ j-1 J \ r-j -I 

To estimate the contribution made by the portion of the sum in the definition of Q r which corresponds 
to the vectors v e R r we obtain by summing over s: 

j-l CO 

pjtw < E p j p u \ R r,jj 



vGR r j s — 1 u—r—j 

r j-l 

< 



j=2 s=l V 7 V 7 u=r-j V 

-- '3-1 \ W„«,-,/r-l\ / P 



J=2 °— 1 V 



3 - 1 / V 1 - P 



The corresponding term for the actual expected values of r\% where v is in the rare set is bounded by 
assumption (II). Hence we obtain 

]T ]T (Ey + ^ (;) )<s r , 

j veR r ,j 



where 



for some c 3 , and 71 = 71 + 5/3, 72 = 72 + 

(iii) Since the sum for Q r contains many terms that cannot be paired with terms in the sum of U r let 
us look at those combinations of r, j, and u that do not correspond to vectors in G r . Let us denote 
by V r .j, u the number of elements for the values of r,j,u that occur in the representation of Q r and are 
not in G r . Those overcounts occur for overlaps which have lengths > M. Denote by iii, 112, ■ ■ • , u r -j 
the individual overlaps of these r — j fictitious 'immediate returns'. Then, if the first intersection is of 
length mi (> M. ) j then 

u-(r-j) 

Vr,j,u ^ ^ ^ ^ \G r J^ u — u r \ , 

M 

Ul = — 

Hence (for u > ^ + (r — j '• — 1) and j < r) 

u—(r—j — l) / _ v / v • , v , M 

u — u\ — 1 \ / r — 1 \ t j ( r — 1 \ 1 u — 



V rJ , u <r- J2 



r- 



J' ! - m \ r-3- 1 J V 3 ~ 1 J jl V J - 1 

where we used the identity J2y=b (b-i) = (b)- Let us P ut ^ f° r tnc set of ficticious vectors v 
that have individual overlaps Ui (i ^ i\, . . .ij) that are not allowed, i.e. where at least one overlap is 
longer than — . Then we can estimate the contribution one gets by counting over H r as follows: 



V(r) = ]T ?P+ 



u(v) 

v£H r 



G 



r— 1 oo 

= E 

r-1 
< ?•' 



J — 1 M=-^+7 J — 1 



r-1 



m 

r - j - 1 



fl j t 3 ( r-l \ m_i / p + 
where we used the identity X)^L a P" ( a - l ) = (i^p) ' Hence 



/r-1 



rp+ Q r (t,p + ). 



fri i! V i- 1 / \i-p+. 

Combining the three estimates (i), (ii) and (iii) yields 

\U r - Q r {t,p)\ < (rp^- 1 + </>) Q r (t,p+) + c 5 q(r - l)Q r _i(i,p+) + S r 

(B) The difference <p(z) = ip(z) — g P {z) between the two generating functions splits into two parts, 
ip = ipi + ip2, which we analyse separately {ip\ reflects the estimates of parts (A-i) and (A-iii), and ip2 
reflects the estimate of (A-ii)): 

(i) The function fi(z) is majorised by the power series 



c 6 (p™ + q + 4>) E \ z ~ 1 \ r rQr(t,p+). 



The sum over r is equal to ^e* 1 p + p +™ (where w = \z — 1|) which can be bounded by 4c* 1 p + 
if \w\ < 1 . A Cauchy estimate with |z| = i? now yields 



£i = 



^(0) 



< g l-2p + -p + R 



If t is large so that f < ^ then we can take R = j and obtain for instance that (assuming P+ < ^) 



using Stirling's formula, where e = c 6 ™ + g + 4^ ■ If £ is small so that for instance f > ^ then we 

take i? = ± and thus obtain £i = ^ </4 r) (0) < c 8 e(2p) r e*^l . 

(ii) The second error function if2(z) is majorised by the power series (t' — j3r) 

00 r J — 1 



= exp 



j=2 s= 





1 N 


)5( 













v 3- 





— exp ■ 



t'i\z- 1| 



1 - (71 +72)71-2 - 1| 1- 727N -1| 

i' 7 k-l| 



< 6771IZ — 1| exp 



1 - (71 +7 2 )7|z- 1|' 



if 1 2 — 1 |7(7i +72) is small enough (e.g. < \), where we have used the identity 



00 r 17-1 / ' 1 \ x s ( r 1 

E E E ( s - 1 ) 1 .7-1 



r=2 j'=2 s=l 



j/- 7 s z r 3 = exp 



1 -y- z 



exp 



1 - z 
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(develop into a Taylor series with variable x and use the identity 5Zfe°=£ y \~ \ y 3 
Hence if we put \z\ = R (R > 1) then 



k—l 



(1 - &)"*). 



Eo = - 



4 r) (o) 



< 



C9771 



3 2 7 iit' 
R r ~ X 



if we assume that (R + 1)771 is small enough (e.g. < |). If i? 
formula) 



then we get (using Stirling's 



1 

r! 



< cxotfiTr 2 ?- 1 - 



(if i' is close enough to t). 

Since P(£ = r) = prV^(O) an d e~*P r (i,p) = <7p (0) we get by combining the estimates (i) and (ii) 



|P(C 



e -'p r (i,p)l < 



^W(O) 



<Ei+E 9 



(and 1 + log 2 < 2) from which follows the result of the proposition. 
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In the following we will apply this proposition to situations that typically arise in dynamical systems. 
There the stationarity condition (I) of the proposition is implied by the invariance of the measure. The 
random variables rjj will be the indicator function of a cylinder set pulled back under the jth iterate 
of the map. Condition (II) is then implied by the mixing property (see below Definition [5]) . The most 
difficult condition to satisfy is (III) because it involves 'short range' interaction over which one has 
little control and which require more delicate estimates (see Lemma [5] below) . A simpler version of 
Proposition!]] is the following corollary (m = 1) which is easily deduced by putting 72 = 0, 7 = l,7i = e. 

Corollary 2 Let M, r be as above. Let rjj, j = 1, . . . ,r, be 0, l-valued random variables and r)jg = Y\ i r\ Vi 
for v € G r . For some S > let R r — Uj=i Rr,j> where R r .j — {v S G r .j : A(w) < S}. Assume \i be a 
probability measure on Q which satisfies the following conditions (Co is a constant, e > 0): 
(L) E(?7j) = (3 for all j = 1, . . . , r (invariance of the measure); 
(LL) Suppose there is a p £ (0, 1) so that for all v G G r j >w \ R r ,j'- 



for all r, j and w; 
(HI) 



|E(Tfr)-pV| <sp w ; 



E (W) < e. 

v£R r 



for all r — 1,2,.... 

Then there exists a constant C\ so that for every t > one has ((, — J2]=i Vj an d t = (1 ~ p) T f3) 



C (C = r) - e-*P r (t,p)| < Ci(e + S0)t r - 1 ^- + C\ (p 




if t> \pr 
(2p) r e*^§ if t < \pr 



3 Measures that are ((f), /)-mixing 

Let T be a map on a space fl and /i a probability measure on Q. Moreover let A be a measurable 
partition of and denote by A n — Vj^o T~^A its n-th join which also is a measurable partition of £1 
for every n > 1. The atoms of A n are called n-cylinders. Let us put A* = U^i -A n for the collection 
of all cylinders in and put \A\ for the length of an n-cylinder A E A* , i.e. \A\ = n if A S A n . 

We shall assume that A is generating, i.e. that the atoms of .4°° are single points in il. 

In the following definition we generalise the 'retarded strong mixing condition' (see e.g. [24]). We 
consider mixing dynamical systems in which the function (f> determines the rate of mixing while the 
separation function f specifies a lower bound for the size of the gap m that is necessary to get the 
mixing property. 



8 



Definition 3 Assume p is a T -invariant probability measure on Q and that there are functions f and 
4> so that: 

(i) f : N — * N (No = N U {0} J is non- decreasing 

(ii) <f> ■ No — * K + is non-increasing. 

We say that the dynamical system (T, /i) is (</>, /)-mixing if 

\H(U n T~ m ~ n V) - fi(U)fi(V)\ < </>(m)fi(U)ti(y) 

for all to > f(n), n > 0, measurable V (in the a-algebra generated by A* ) and U which are unions of 
n-cylinders. 

Systems that are (</>, /)-mixing are for instance: 

1. Classical ^-mixing systems (see, e.g. [5]): / = 0. These include equilibrium states for Holder 
continuous potentials on Axiom A systems (which include subshifts of finite type) or on the Julia 
set of hyperbolic rational maps. In this case the partition A is finite. 

2. Dispersing billiards [23]: / is linear. 

3. Equilibrium states for Holder continuous potentials (that satisfy the supremum gap (see section 5) 
on the Julia set of rational maps where the Julia set contains critical points: / is linear, <f> is 
exponential. 

4. Multidimensional piecewise continuous maps [22]: / depends on the individual cylinders (|„4| < 
oo). 

For r > 1 and (large) r G N let as above G r (r) be the r- vectors v — (vi, . . . , v r ) € Z r for which 
1 < vi < V2 < ■ ■ ■ < v r < t. Let t be a positive parameter, W C 51 and put t = [t/fi(W)] be the 
normalised time. Then the entries Vj of the vector v € G r (r) are the times at which all the points in 
CV = rfj=i T~ Vj W hit the set W during the time interval [1, r]. 

Lemma 4 Let (T,/i) be ((f), f)-mixing. 

Then for all r > 1, Wi C f2 unions of rii- cylinders, i = 1, . . . , r (n, > 1), and all 'hitting vectors' 
v G G r (r) with return times u^+i — Vi > f{n,i) + rij (i = 1, . . . , r — 1) one has 



{phi .mixing} 



n; =1 urn 

and d(v, ft) = van\i{vi + i — m — ni). 



< (l + <l>(d(v,nW-l, 



A consequence of this is that there exists a < r/ < 1 so that for all fi(A) < ry ' for all A 6 A* . 
3.1 Estimate of the rare set 

In this section we provide an estimate for the rare set for general ((f), /)-mixing maps. We will then 
use this result in its full strength later to show that the return times distribution at periodic points is 
compound Poissonian for rational maps that have critical points. For a 'hitting vector' v € G r (r) (r a 
large integer) we put C$ = Dfe=i T~ Vk W. Let 6 > f(\ W\) (W a union of cylinders of the same lengths) 
then 

Rv,j(t) = {ve G rj (r) : min(v 4fc+ i - v ih - \ W\) < S}, 

k 

where the values v ^ , . . . , Vi . are the beginnings of the j blocks of immediate returns (notation as in 
section O(H)). 

Lemma 5 Assume (T, fi) is (<j), f) -mixing and assume that there is an m € N so that for every n for 
which f(n) < 5 there exists an M < n so that A n n T~ £ A n ^ for £ < M implies that £ is a multiple 

of TO. 



{product .mixing} 



{section. return. 



{R. small} 
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Then there exists a constant C'2 so that for all n-cylinders A r , 

r 3-1 



veR r j=2 s=l ^ ' ' \ " J 

where: 

(i) n',m' (ml < n') satisfy f(n') < M — n' and f{m') < to — ml , 

(ii) 7 > 1 + (f>(mim(vi+i - Vi) - n'), 
(lii) A n , e A n , A n c A n >, 

(iv) A m , e A m ', A n , c A m ,. 

Proof. As in section (A-ii) of the proof of Proposition [TJ put Rf, ^ for those v £ R r j for which 
Vi + i — Vi > S for s — 1 indices i%, . . . , i s -i (i s — w r , s < j — 1). We consider two separate cases: (I) 
s > 2 and (II) s = 1. 

(I) Assume s > 2 and i\, 22, ... , i s -i be the indices for which Vi k +i — Vi k > 5 > f{n) for fc = 1, . . . , s— 1. 
All the other differences are > M and smaller than S. Let A n > be an n'-cylinder so that A n C A„< 
where n' is so that f(n') < M — n' . Let j be the number of blocks (i.e. v £ G rj ) and let i l7 . . . , be 
the beginnings of the 'blocks of immediate returns' (clearly s < j — 1). There are r — j immediate short 
returns of lengths e [m,n). Let us put 

W lk = A n for k = l,...,s, 
Wi> k = A n , for k=l,...,j, 

Wi = A m , n T- m A m , n T- 2 "M m - n ■ • • n T-^-^ m A m , for all t g {i fe : jfc} U {i' fc : A:} 

where «i is the overlap for the ith return (which is an immediate periodic return). By our choice of n' 
we have achieved that 1^+1 — Vi k > 5 > f(n) and Vi+i — Vi > f(n') for i ^ ik, k = X, . . . , s — 1 and 
i ^= i' kl k = 1, . . . , j. By Lemma 0] we obtain 

(r \ r 

p| T-^Wi < a^" 1 H^Wi) < arV(^) J "X^) s (a2M(AnO) u , 
i=l / i=l 

(ai = 1 + min(0(<5 — n),(j)(n — n')), 0:2 = 1 + 0(m — to')) where the components of n = (n\, . . . ,n r ) 
are given by ni k = n for k = 1, . . . , s and rt^ = n' for i ^ fe = 1, . . . , s), where u = Mi is the total 
overlap. 

The cardinality of i?* ■ u — R^, n G r .j :U has been estimated in part (A-ii) of Proposition Q] to be 

rc^('.:i)^(;:i)(r-J±. 

Therefore 

E KCe) < oT 1 ( jl ; ) ^^^(A„0r s (;:;)( p ^ , ) (a 2 M(A ro ,))" 

(II) If s = 1 then all returns between blocks are less than S for all k. In the same way as above we 
obtain 

E moo ^ «rV(40(M4»')) i - 1 f !Jl ! K r "7 \ ) M^mor 



i - 1 / V r _ - 1 



Summing over s and using the estimates from (I) and (II) yields 

3—1 00 

E = EE E E mso 

veR r j s =l 11=1 — j v£R s ■ 



10 



j=2 s=l V / ■ «=r-j \ J / 



j=2 8=1 



The lemma now follows since piaip^fH") ~ a '/ i (^m / ) with a a' which is slightly larger than a 2 - Now 
we write a 2 instead of a'. 1 

In the case of classical ^-mixing maps (see subsection 3.2 below), when / is zero, we get the following 
simpler result. (We simply put nl = n and m' = m which then results in A n > — A n and A m / = A m .) 

{R . small . phi-mix 

Corollary 6 Assume (T, /x) is (^-mixing and assume that there is an m € N so that for every n there 
exists an M < n so that A n n T A n for I < M implies that £ is a multiple of m. 
Then there exists a constant C3 so that for all n-cylinders A n : 



veR T j=2 s=l 

where a = 1 + 0(0) and A m G A m contains A n . 
3.2 0-mixing measures 

We say that the dynamical system (T, fx) is 4>-mixing if / is identically zero, i.e. 

| M (t7 n T- m -"F) - /x(f/)/x(y)| < 0(m)n(U)n(V) 

for all m, measurable V (in the cr-algebra generated by A*) and [/ which are unions of cylinders of the 
same length n, for all n. The function 4> is assumed to be monotonically decreasing to zero. 

Let W be a set in 17. Then the entries Vi of the vector v G G r (r) are the times at which all the 
points in C$ = fX = iT~ Vi W hit the set W during the time interval [l,r]. Following Lemma 2] we get 
that for n^-cylinders W% Cfi,t = l,...,r; 



< (1 + ^(d(w))) r - 1, (2) {phi .mixing} 



for all 'hitting vectors' v G G r (r) with return times v i+ i — > n; (i = 1, . . .r — 1) where d(v) 
minj(w i+1 - Vi - nA. 

3.3 Distribution near periodic points for 0-mixing measures 
Lemma 7 Let x be a periodic point with minimal period m. If fi is <fi-mixing then the limit 



p = lim 



j\ogn(A em (x)) 



exists. 



Proof. We show that the quantity inside the logarithms is nearly superadditive. Let A be an integer 
so that 4>(Am) < h. Then we have 

\\o gf i(A km+Am+em (x))\ > \\o gf i(A km nT- km - Am A em (T- km - Am x))\ 

and by the mixing property 

fi(A km n T- km - Am A hn (x)) = n(A km )fi(A lm (x))(l + O'^Am))) 
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(where O* means that 



< 1 for all e). If we put aj = \log n(Aj m (x))\, then 

ak+A+e >a k +a e - |log(l - <f>(Am))\ > a k + a £ - 20(Am) 
for all positive integers k, £. Iterating this inequality yields 

«rfc+(r-i)A+s . ra k - 2(r - l)^(Am) 



(rfc + (r — 1)A + s)m 



> 



> 



(rk + (r — 1)A + s)m 

1 a k 2<p(Am) 



1 



A 

k 



kr 



km 



for positive integers k, r and s £ [0, fc + A). If we put n = kr + (r — 1)A 
let r — ► oo we obtain 

a„. . 1 a k 2<j){Am) 



A ' 

- s, < s < k 



A - 1, and 



lim inf — — > 

n— >oo nm 1 



Now let k —> oo and we finally get 



which implies the lemma. 



lim inf > lim sup - 

n— >oo nm k^oo km 



k + A ' 



(ik 



I 



As a consequence of the lemma we see that p < r\ m for some rj < 1. (This follows from the fact that 
m-cylinders have measure < rj m for some T) < 1 |13j.) In particular p is always strictly less than 1. 

In the following we shall assume the stronger property that p = linin^oo ~^x7-ifr • This of course 
implies the limit in the lemma, but we are not sure whether the reverse implication is generally true. 



Also put q n = swp i>n 



fj,(A i + m (x)) 



For t > and integers n we put Q f°r the counting function 



Sj=o XA n (x)) ° T J with the observation time 



(where x is periodic with minimal period m) 



(l-p)n(A n (x)) 



In order to satisfy the assumptions of Proposition [T] we put 7 = a, 71 = a5 n jj,(A n ) and 72 = a/u(A m )) 
and Corollary [SI 

Theorem 8 Lei (/i, fl) be a <j)-mixing measure with partition A (finite or infinite), x a periodic point 
with minimal period m and p and q n as above. 

Then there exists a constant C4 so that for every S > and every t > one has 



„2r 



t r e 2r+§t 



if t > \pr 



|P(C* = r) - e-*P r | < (7 4 ^(A n )r- 1 — + C 4 (p» + g„ + { J"' > r t i±|E :, , , - 

n I (2p) e 1 ~ 4 *' «/ t< T^pr 

Proof. We use Proposition [T] and have to verify conditions (I)— (III) . From the definition of p and q n 
assumption (I) is clearly satisfied with p± = p ± q n . 

To verify condition (II) let v G G r j. u and let us look at the measure of /i(CV). By ([2|) we have (A 
is as defined in section [ 



{phi-mixing} 



tic*) - n 



k=l 



<((l + 4>(A(v)-n)y -l)]Jv(D k ), 



where D k is the fcth block, i.e. 



D k = p| T-^A n {x). 



k=l 
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Since fi(Dk) = p(A n+mUk (T Vz >* (x))) we get by definition of q n 

lAPk) K A n+mu k ) 
H(A n ) K A n) 

) n(A n+mUk ) 



and therefore 



p,(A n ) fi(A n+m ) ^(A„ +m ( Ufc _!)) 
(p + 0(q n )r 



ri !#y = n (p + °^))" fe = (p + • 

fe=l ^ k=l 



where u = u k- Hence 



.7 



< At (^(( P + g n )«-^ 



fe=l 

and consequently 

| M (CV) -pX^'l < n(A n y ((p + g n ) u - P u + 0(A(tf) - n)) 3 - 1)) . 

Hence, if v £ R r j then we get assumption (II) with 7 = a, p± — p± q n (and 7(71 + 72) < if m, ; 
are not too small). Here we use M = n — m. 

To verify assumption (III) we use Corollary [6] . We obtain 



where a = 1 + 0(0). Hence condition (III) of Proposition [T] is satisfied with 71 = 8p,(A n {x)) 1 72 = 
afj,(A m (x)) and (3 = n(A n (x)). I 

Let us note that this result applies to finite as well as infinite partitions A. Since here we focus on the 
recurrence properties around periodic points we do not require the condition J2aeA ^°&^{A) < 

00 (which is necessary in order to get finite entropy or the theorem of Shannon-McMillan-Breiman) . 

Equilibrium states for Axiom A systems: Let us now assume that p is an equilibrium state 
for a Holder continuous function /(with pressure zero) on an Axiom A space (shift space) which has 
the finite, generating partition A (see [5]). Then p = hv where h is a normalised eigenfunction for 
the largest eigenvalue of the transfer operator and v is the associated eigenfunction. In particular v 
is e~^-conformal, i.e. if T is one-to-one on a set A then v(TA) — j A e~f dv{x). If we replace / by 
/ = / + log h — log h o T then p is e~^-conformal. Thus, if x is a periodic point with period m, then 

p(A n (x)) = p(T m A n+m (x)) = f e-? m M d p(y) = /i(A n+m (x))g n e-™, 

where q n is a number that can be estimated by 

I log Qn I < var„/ m < const. (var„/ + var n+m log h + var n log h) < const. 0™, 
for some 9 € (0, 1) (Holder exponent). Hence 

p{A n+m (x)) _ 

fl(An(x)) 

where p = ^ and q n = p(q n — 1) can be estimated by |g„| < p9 n const.. In particular the limit 
linin^oo exists and equals p. It is known that p is 0-mixing where <p(k) = p k for some 

p € (9, 1). Let us now apply Proposition 2 and in order to minimise the term a n — C\ ip™ -t- q n -\- p 11 ) 
we choose S„ — l °sK A n(. x )) _ Then e n < const. (p» + np(A n (x))) (again M = n — m). 
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Corollary 9 Let fi be an equilibrium state for a Holder continuous function on an Axiom A system. 
Then there exists a constant C5 so that for all periodic points x, t > and r = 0, 1, . . . one has (p is as 
above): 

|P(C* =r) -e"*P r | < C^A^llog^A^xW-^+Cs (p* + n»{A n {x))) 



(2 P y 



2r+%t 



g 1 — 4p 



i > 
t < 



\pr 
\pr 



Algebraically </>-mixing systems. If we assume that /1 is 0-mixing (with respect to the partition A) 
where <fi(k) — (D(k~ K ) for some k > 0, then let us note that 

n{A n+m {x)) {1 + 0(O))Li(A n (x))fi{A m (x)) 

p + qn= ,(A n (x)) ~ KMxj) " ClKAm(x)) 

implies the very rough estimate q n < /j,(A m (x)) . With S n — [i(A m (x))~ « one now obtains (n >> m) 

Cn < c 2 (p» + f*(4 m (i)) + o- K ) < c 3 ^(A m (x)). 

Corollary 10 Let /j is cf>-mixing and 4>{k) ~ /or some k > 0. TTien i/iere exists a constant Cq so 
that for all periodic points x, t > and r = 0, 1, . . . one /ias is as above): 

* = r) - e-*P r | < C 6M (A n )| log/i^Cx))^" 1 ^ + C 6 n M (A m (x)) { ^^ + t ^ t ! J 1^ ■ 

r\ I (2p) e !- 4 p i/ t < T^pr 



3.4 Example 

In [T9l [20] it has been shown that for ergodic systems every possible distribution can be realised for 
entry and return times of ergodic systems if the sequence of sets is suitably chosen. Naturally all 
settings in which the limiting distributions are shown to be exponential or Poissonian (in the case of 
higher returns) have to assume that the target set is a cylinder set (or a topological ball as in [231 fTTj ). 
Here we show that even if we take cylinder sets then there are points which do not have a limiting 
distribution at all. 

For simplicity's sake let S be the full two shift with symbols 0, 1 on which we put the Bernoulli 
measure with weights w, 1 — w > (w 7^ |). Let y — 0°° and z = 1°° be the two fixed points under 
the shift transformation a. They have periods mi = m-i — 1. The entry times at y,z are compound 
Poissonian with the p- weights p\ = w and pi = 1 — w. Put e = -^\px — pz\ and we will now produce a 
point x so that the return times distribution up to some order ro oscillates between the two compound 
Poisson distributions. Choose n\ so that the cylinder A ni (y) = A ni (0 ni ) has the distribution 

|P(C= r )-e-*Pr(t,Pi)| < I 

for t < to and r — l,...,r n for some to > 0. Now we choose n 2 > n i so that for the cylinder 
A n2 (0 n A n ^ ni ) one has 

|P(C* 2 =r-)-e-*P r (*,p 2 )| < I 

for t < to and r = 1, . . . ,ro- This can be done because the limiting distribution is invariant under 
the shift a (i.e. the limiting distribution of the cylinder A n2 (0 ni l™ 2_ni ) as ni — * 00 is equal to the 
limiting distribution of the cylinder A n (l°°) as n — > 00). Continuing in this way we find a sequence 
of integers n\, ri2, n^, . . . so that the distribution of alternates within an error of | between the 
distribution e~ t P r (t,pi) (for odd j) and e~ t P r {t,p2) (for even j) for t < t and r < ro- Hence the 
point x = f)^ A nj (o ni l n2_Tl1 • • • ™i) (* is is j is odd and 1 if j is even) has no limiting 

distribution. 

Naturally, this construction can be carried out in all </>-mixing systems. Instead of two fixed points 
one can also take any finite number of periodic points and then construct a point which takes turns 
visiting all of those so that at each visit it stays long enough so that its return time distribution gets 
arbitrarily close to the return time distribution of the periodic orbit it visits. 
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4 Return times 

Instead of looking at the probability of a randomly chosen point in the space f2 to enter a given set A, 
here we look at the statistics with which points within A return to A again. In the case of the first 
entry and return times, these two distributions have for general ergodic systems been linked in |12j . 
Higher order entry and return times have been related in [4]. It turns out that these distributions are 
the same only if the first return time is exponential. Similarly, the number of entry and return times 
have the same distribution if it is Poissonian. However, near periodic orbits we get for the return times 
a distribution which is very similar, namely it is in the limit given by the following compound Poisson 
distribution. 

Let p £ (0, 1). If we define 



{section. return. 



P r (<,p)=]>>'- J (l 

3=0 



P) 



for r = 1,2,... and Pq = 1 — p then the generating function for the probabilities e t P r is 



r=0 



zp 



The mean of this distribution is |^ and the variance is Tfz^pr ■ Again note that if p = then we get 

the Poisson terms e~ t P r (t, 0) = e~ 1 ^ and the generating function e* 1 - 2-1 *' which is analytic in the entire 
plane whereas for p > the generating function g p (z) has an essential singularity at -. The expansion 

at zq = 1 yields g p (z) = Y^=q{ z ~ l) fc Qfc where 



)k(t,p) 



{l-pfj- 



1=0 



(Qo = 1) are the factorial moments. 

For a set A let us now define the random variable C,a 



Xa Ej=i XA ° T j and put Q = C,A n (x) where 



t = (1 — p)r n /x(j4 n (a;)); we also denote with fi n the conditional measure to the cylinder A n (x). In a 
similar way we can now prove the following result. 

Theorem 11 Let (/i, f2) be a (j)-mixing measure with partition A, x a periodic point with period m and 
p and q n as above. 

Then there exists a constant CV so that for every S > and every t > one has 

2r ( t r 2r+— t 'f -f \ 1 

P(C* =r\A n )-e- t P r <C 7 n5^A n )t r - ie — + C 7 n(p^+q n + m)\ it* Z ZT ' 

r! I (2p) r e 1 - i p if t < ±pr 



{phi-mixing} 



where r„ = 



(l-p)p,(A n (x))- 



If we compare these error terms to the ones for the entry times, we notice the additional factor n which 
comes from satisfying the condition (I) of Proposition [1] (cf. [T3]). 

Let us note that for r = this result has previously been obtained by Hirata [15] for equilibrium 
states for Holder continuous function on Axiom A systems. Here however we also get error estimates: 

P(C* = 0|4») ~ (1 -P)e~* < (p™ + nn{M*))) ■ 



Note that if p > then Pq(Q,p) = 1 — p is strictly less than one and P r (0,p) = p r (l — p) for r > 1. 
There is a point mass at t — which corresponds to immediate returns within the neighborhood of the 
periodic point. These are clearly geometrically distributed. 
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Remark: By adapting a recent remark of Chamoitre and Kupsa [4], we proved in [14] under the 
condition of the existence of the asymptotic distribution of successive return times that the asymptotic 
distributions for the entry and return times are related by the formula (k = 1, . . .) 

D k (t) = J (p k -i(s) - Dk(xj) ds 

where Dk(t) is the limiting distribution P(££ = k) as n — > oo, and Df.(t) = ]im n -> 00 P(Q = k). 



5 Rational Maps 



Let T be a rational map of degree at least 2 and J its Julia set. Assume that we executed appropriate 
branch cuts on the Riemann sphere so that we can define univalent inverse branches S n of T™ on J for 
all n > 1. Put A n = {f(J) '■ <f £ S n } (n-cylinders). Note that the diameters of the elements in A n 
go to zero as n — > oo. Moreover, A n is not the join of a partition, yet they have all the properties we 
require. 

Let / be a Holder continuous function on J so that P(f) > sup / (P(f) is the pressure of /), let fi be 
its unique equilibrium state on J and £„ — 2j=i XA n ° P^ the 'counting function' which measures the 
number of times a given point returns to the n-cylinder A n within the normalised time t„ = [t/ fi(A n )]. 
Although n is not a Gibbs measure we showed in [13] that for almost every x 



P(C« - r) 



r 



as n 



Theorem 12 Let T be a rational map of degree > 2 and \x an equilibrium, state for Holder continuous 

f (with p(f)> sup /;. 

Then there exists a p £ (0, 1) and Cs so that for every periodic point x £ J the return times are 
approximately compound Poissonian with the following error terms: 



\nC = r) 



e-*P r \ <C & ~p n t r - ie ^- 



C 8 ^ 



if t>h 



pr 



(2tr) r ev i- 2 ? if t < \pr 



wh 



ere p - 



e ' ^ mP W) and m is the minimal period of x. 



The univalent inverse branches S n of T n (with appropriate branch cuts) split into two categories, 
namely the uniformly exponentially contracting inverse branches S' n and the remaining 5" = S n \ S' n 
for which do not contract uniformly. In |llj we showed the following result: 

Lemma 13 Lemma 9) Let r\ £ (0, 1). Then there exists a constant v > so that for all r > 1 

and v — (vx,V2, ■ ■ ■ ,v r ) £ G r satisfying mirij(vj +1 — Vj) > (1 + v)n (clearly r < iA j: 



- i 



< 



for all sets W\, . . . , W r each of which is a union of atoms in A n and for all large enough n. 

Let us define the rare set R r : We put R r for the set all v £ G r (r„) for which minj(vj+i—Vj) < (l+q)n. 

Lemma 14 Let x £ J be a periodic point with (minimal) period m. For all large enough n one has 
that A n (x) n T~ l A n (x) ^ for I < n/2 only if £ is a multiple of m. 

Proof. Put n = km + n', < n' < n, and tp — ijj k ■ ■ -ip 1 ^ 11 , where ip 1 , ...,tp k £ S m , <p n £ S n >. 
Since x £ A is periodic with period m we get that T lm A n A ^ and in particular x £ T" n A for all 
i = 1, . . . , k. Since the sets ip(J n O m ) are all disjoint for different i, we obtain ip 1 = ip 1 for all i. Put 
ip — ip 1 and we get <p = (ip) k cf> k (with ip concatenated k times). 



{rational .maps} 
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{rational . period 
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Now assume that A n T~ l A ^ for some i < \ which is not a multiple of m. Since for some i, 
im < £ < (i + l)m and T vm A n T~ e+tm A ^ 0, we can assume that £ < m. Suppose that there are 
arbitrarily large n so that £ < m and V — A n T~ l A ^ 0. Similarly as above we put n = k'l + n" 
(0 < n" < £) and obtain that S S n decomposes as <f) = (tp) k <p n where V € St, 4> n £ 

Now since (^>) fc (J n fl^) — > x as fc' — > oo, and x is periodic with period m, we see that such £ < m 
cannot exist. Hence, for all n large enough T e A n ^4 ^ and £ < § implies that ^ is a multiple of the 

period m. 1 

Proof of Theorem 1121 We are going to verify the conditions of Proposition [TJ Let x £ J be periodic 
with minimal period m. Then 

(I) holds by invariance of the measure = /J.(A V ) for all j. 

(II) Since /i = hv where ft, is a Holder continuous density and v is e~^-conformal we obtain as before 
that 

/i(4»(a:)) = MT m A n+m (x)) - / e~t m to dfi(y) = ti(A n+m (x))q n e-? m W , 

JA„ +m (x) 

where we have used that fact that fi is e~-' r -conformal with respect to the function / = /+log h— log hoT. 
The factor q n satisfies | logg n | < var n / m < const. 6 n , for some 9 £ (0, 1). Hence 

n(A n+m (x)) _ 

fi(An(x)) 

where p = (independent of n) and the error term g„ = p(g„ — 1) is bounded as |g„| < cip9 n for 

a constant ci which is independent of the periodic point x. 

(III) Here we use Lemma[SJ By Lemma mi we can choose M = [n/2]. Furthermore we set S = (1 + v)n. 
According to Lemma Hcfl our separation function / is given by f(k) = (1 + v)k. Hence n' = [n/ (1 + v)] 
and m! = [m/(l + v)]. Then A n / is the n'-cylinder that contains A n = A n (x) and whose measure is 
fi>(A n i) < ( o"/( 1 + 1 '). Similarly A m / is the m'-cylinder that contains A m (x) and and whose measure is 

n(A m i) < p m /( 1+v \ Let us choose p < 1 so that p > max , rj, ^ . Then (for all large enough n) 

E mcw < 1 g ( i : ; ) tj-^ ( ; : ; ) , 

iTefl r j=2 s=l v ' ■ \ / 

where 71 = 5/i(A„/) < p", 72 < a/i(A m /) < p m , [3 — rp(A n ) and by Lemma [T3l a = 1 + i] S ~ n . 
Moreover, since p + = p + q n < p m , — < p™, f>+ — P- < q n < p n and = 0(<5 — n') < < p n one has 

p+ + p+ - P- + <t> < ap™ 

for some c\. The theorem now follows from Proposition [T] I 
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